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FINITE  WIDTH  CURRENTS,  NON-UNIFORM  MAGNETIC  SHEAR 
AND  THE  CURRENT  DRIVEN  ION  CYCLOTRON  INSTABILITY 


In  our  recent  study^  on  the  effects  of  magnetic  shear  on  the  current 
driven  ion  cyclotron  instability  (CDICI),  we  obtained  the  intriguing  result 
that  even  a  small  shear  could  produce  a  significant  reduction  in  the  growth 
rate.  The  transition  from  the  local  to  the  non-local  treatment  essentially 
induces  a  singular  perturbation,  (the  small  shear  being  the  coefficient  of 
the  highest  —  here  the  second  —  derivative  term),  and  there  is  also  a 
change  in  the  boundary  conditions,  from  plane  waves  to  a  bounded 
solution.  It  was  thus  possible  to  understand  and  justify  the  sharp  change 
in  the  growth  rate  from  the  local  to  the  non-local  treatment. 

However,  it  was  assumed  in  the  treatment  above  that  both  the  shear  and 
the  current  were  uniform  in  space.  In  laboratory  as  well  as  space  plasmas 
on  the  other  hand,  one  necessarily  deals  with  finite  width  currents,  so 
that  the  scale  length  of  the  current  profile  is  an  additional  parameter 
which  could  play  a  rule  in  determining  the  precise  nature  of  the  transition 
from  the  local  to  the  non-local  treatment.  One  can  see,  a  priori,  that  if 
the  current  width  scale  length  Lc  is  much  larger  than  the  shear  length  Ls, 
the  effects  due  to  shear  should  prevail  unhindered.  In  the  other  limit, 
when  Lc  is  much  smaller  than  Ls>  it  is  also  clear  that  the  total  change  in 
the  angle  of  the  magnetic  field  vector  in  traversing  the  slab  will  be  so 
small  that  the  non-local  effects  due  to  shear  may  not  come  into  play,  and 
the  local  growth  rate  should  prevail,  apart  from  any  non-local  effects  due 
to  the  variation  of  the  current  profile  itself.  Thus  a  transition  from  the 
local  to  the  non-local  results  will  occur  as  the  current  width  is  increased 
from  a  value  much  smaller  than  the  shear  length  to  a  value  larger  than  the 
shear  length;  this  transition  can  be  expected  to  be  primarily  governed  by 
the  ratio  of  these  scale  lengths,  and  a  quantitative  study  of  this 
phenomenon  is  the  main  objective  of  this  memo  report. 

An  additional  feature  to  be  considered  is  the  relationship  between  the 
current  profile  and  the  shear.  If  the  shear  is  entirely  generated  by  the 
driving  current,  we  will  call  it  the  self-consistent  shear.  If  the  shear 
is  generated  by  (stronger)  external  currents  which  do  not  drive  the 
instability,  we  will  call  it  the  prescribed  shear.  These  two  cases  could 
Lead  to  slightly  different  results.  Examples  of  one  situation  or  the  other 
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can  be  found  in  various  plasmas  of  practical  importance,  and  there  are  also 
situations  where  the  total  shear  derives  from  both  the  mechanisms • 

With  the  introduction  of  a  finite  slab,  the  problem  can  no  longer  be 
translated  to  a  convenient  origin  where  the  given  wave  vector  is  perpen¬ 
dicular  co  the  local  field  line.  Instead,  we  will  take  the  center  of  the 
slab  as  the  origin  and  introduce  explicitly  k°,  the  parallel  component  of 
the  wave  vector,  as  a  prescribed  parameter.  One  must  maximize  over  all 
possible  k°,  even  in  the  non-local  theory,  to  obtain  the  optimal  normal 
mode  solution.  We  introduce  a  finite  current  width  along  the  x  direction 
by  taking  the  electron  distribution  function  to  be  a  shifted  Maxwellian 
with  a  drift  velocity  parallel  to  the  magnetic  field  to  be 

-2  v 

V  (x)  =  V°  g(x  /L  ),  g(C)  =»  e"^  ,  -  x  +  ^  .  (la) 

a  age  g 

The  Vlasov  equation  is  satisfied  since  xg>  guiding  centre  position  is  a 

constant  of  the  motion,  v^  being  the  electron  velocity  component 
peroendicular  to  the  magnetic  fiaLd  direction  z  and  £  is  the  electron 
gyrof requency .  A  very  small  correction  in  the  constant  of  motion  due  to 
shear  can  be  ignored  here.  The  current  profile  is  then  found  to  be 


j(x)  «  n  eV°axp[-x2/(L2  +  p2)]  *  n  eV°g(x/L  ), 

rt  n  *  a  '  O  Q  £ 


where  «  L2  can  be  ignored;  this  is  also  equivalent  to  replacing  Xg  by  x 

in  equation  (la).  This  will  generate  a  self-consistent  shear  in  the 

magnetic  field,  given  in  terms  of  5  *  —  by 

c 

,  x  4TTn  e  V°L  5 

B  (x)  -  —  /  j(x)dx  *  - - -  /  g(C)dC,  (2) 

y  o  o 

B  (x)/B  -  (L  /L  )  f  3(S)dC, 

7  2  c  s  * 


(1/s)  -  L 


Atth  eV 
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and  a  corresponding  variation  of  the  parallel  wavenumber, 


k  (x)  -  k°  +  k  3  (x)/B  , 
t  J  7  7  2 


2 


or 


u(x)  =  k  (x)/k  »  u  +  (L  /L  )  /’  g(5)d€.  (4b) 

II  y  o  c  s 

J  o 

The  shear  length  is  determined  by  the  physical  parameters  in  eq.  (3). 

Following,  the  usual  procedure  for  introducing  the  non-local  effects, 1 
we  let  +  -(l/i)d/dx,  and  assuming  that  the  higher  derivative  terms  are 
less  important,  we  obtain 

[p?  +  Q(u(x),  V  (x)]j(j>  =  0,  (5) 

1  dx2  d 


or 


where  o , 

Of  ki;,  U 
by 

Q  *  -  P^Qj/A  (&) 

where  0^  and  A  are  described  by  eqs.  (6)  and  (7)  of  Ref*  1.  In  general,  it 
would  be  necessary  to  solve  (5)  by  a  numerical  method  such  as  the  Numerov 
shooting  code  technique.  The  case  of  pure  shear,  without  any  effects  due 
to  finite  current  width  (V^(x)  =  V°),  was  discussed  in  Ref.  1  by  expanding 
Q  to  second  order  around  uQ>  the  angle  of  maximum  growth  in  local  theory 
given  by  Q(u0)  m  0*  The  resulting  Weber  equation  provides  a  dispersion 
relation  for  the  determination  of  the  eigennode  frequency.  When  anharmonic 
terms  become  important,  the  error  involved  in  approximating  Q  by  a 
quadratic  form  around  uQ  will  become  significant.  An  improved  procedure  is 
to  consider  an  expansion  around  the  (complex)  position  where  the 
derivative  of  Q  with  respect  to  5  vanishes, 

q'(w,  5l)  -  0.  (7) 


:iir>2  Vi*  ■  °- 

c  d 5 


is  the  ion  Larmor  radius.  Apart  from  the  change  in  definitions 
and  Vdf  to  include  the  space  variation  of  the  current,  Q  is  given 
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Then  the  dispersion  relation  (analogous  to  eq.  (18)  of  Ref.  1),  reduces  to 
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Q(,jJ >  CL)  -  (24  +  1)  (Pi/Lc]  [-  -j-  Q”(<^»  f.^]172,  (8) 

where  i  is  the  mode  number.  Equations  (7)  and  (8)  are  two  complex 
equations  for  the  simultaneous  determination  of  the  complex  eigenf requency 
a)  and  the  complex  angle  u^,  corresponding  to  the  complex  position  % ^ .  We 
have  tested  the  results  obtained  from  eqs.  (7)  and  (8)  against  the 
numerical  shooting  code  solution  of  eq.  (5)  in  a  variety  of  cases,  and  the 
results  are  in  excellent  agreement  for  the  eigenf requencies  as  well  as  the 
solutions  for  the  potential  <f>. 

Considering  Q(u,  v)  as  a  function  of  two  variables  u  and  v,  each  a 
function  of  the  explicit  expressions  for  the  derivatives  are 


<r  ■  V*v«  %  -  K'.ty'K' 


Q"  -  'L  /L  ; 2 s2 Q  -  2'L  /L  SgQ 

c  SJ  UU  "  C  S'  u 


-  4  (l  /l  }(v°/v  ) eg2 Q 
^  c  d  e;  uv 


-  2g(l-252)(V^/ve)Qv  +  4(V°/vJ2r2g2Q 


d'  e' 


vv 


(9a) 

(9b) 


We  have  solved  eqs.  (7)  and  (8)  for  a  range  of  slab  widths  (L^/p^)  3 

LO  to  107  and  shear  lengths  (L  /p  )  *  10^  to  10^.  The  results  are 

s  1 

displayed  for  l  *  0  mode  in  Figures  1-5.  The  basic  transition  from  the 

local  to  the  non-local  results  is  evident  in  Figure  1,  which  depicts  the 

growth  rate  in  units  of  ion  Larmor  frequency,  (y/^)  as  a  function  of  the 

width  of  the  current  slab  expressed  in  units  of  ion  Larmor  radius, 

(L^/p^).  Typical  parameters  are  the  ion  to  electron  mass  ratio 

u  *  1837,  the  electron  drift  to  electron  thermal  velocity  ratio 

V^/v^  »  0.28,  the  ion  to  electron  temperature  ratio  t  »  (T^/Te)  =*  0.5  and 

the  transverse  wave  number  given  by  b  =  (1/2)  -  0.6.  Three  different 

7  i 

shear  lengths  are  represented  by  (L  /p  )  *  10^,  1 0 5  and  10^.  In  each  case, 

when  Lc  i3  sufficiently  large,  one  obtains  the  non-local  growth  rate  y/f.^  * 
0.0127.  Since  the  shear  strengths  considered  here  are  all  small,  we  find 

the  same  non-local  value  for  y.  For  a  stronger  shear  (p  / L  )  *  10~2,  (not 

1  s 
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Y 

Figure  1  A  plot  of  the  growth  rate  (— )  against  the  current  channel 

width  (L  /p  )  for  three  different  magnetic  shear  lengths . 
c  i 

Here  *  0*28  v  ,  p  *  1837,  b  ■  0.6,  and  t  *  0.5. 


01  0.001  0.01  0.1  1.0 

(L-c/Ls)— - 

A  plot  similar  to  Figure  (2a)  except  that  the  value  of  b  here 
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A  plot  for  the  complex  angle  for  various  values  of 

Lc/Lg.  Here  b  =  0.6,  y  =  1837,  V°  =  0.28  v~,  and  t  *  0.5. 


A  plot  of  the  growth  rates  against  b. 
Here  T  =0.5,  u  =1837,  V°  =  0.28  ve*  The  top  solid  line  is  a 
solution  of  the  local  dispersion  relation.  The  second  solid 
line  is  for  Lc  =  =  106p^t  The  dotted  lines  show  the 
transition  from  local  to  nonlocal  values  as  Lc  approaches 
Ls.  The  lowest  solid  line  is  for  L„  =  L„  *  10^o . . 


shown  in  Fig,  1),  a  lower  non-local  value  is  obtained,  (y/ft^)  =  0.0115,  and 
still  lower  non-local  growth  rates  would  be  attained  if  the  shear  is 
further  increased.  In  the  other  limit,  when  Lc  is  sufficiently  small,  the 
local  growth  rate  prevails,  (y/2^)  =*  0.0296.  The  transition  curves  for 
different  shear  values  look  very  similar,  apart  from  a  shift  in  Lc>  and 
indeed  as  shown  in  the  next  figure,  possess  a  universal  behavior  as  a 
function  of  Lc/Ls.  Another  interesting  feature  in  Fig.  1  is  the  decrease 

n 

in  growth  rate  for  Lc  <  10^p^  which  arises  from  the  nonlocal  effect  of  the 
space  variation  of  the  current  itself.  The  slab  width  is  too  small  to 
obtain  any  significant  variation  of  the  direction  of  the  magnetic  field  due 
to  shear,  and  thus  this  reduction  in  growth  rate  is  independent  of  the 
shear  parameter.  If  the  current  width  is  further  reduced  to  just  a  few  ion 
Larmor  radii,  the  instability  is  completely  quenched.  We  can  call  this  the 
filaments!  quenching  —  current  filaments  only  a  few  ion  Larmor  radii  in 
width  cannot  sustain  the  instability.  This  feature  may  be  of  practical 
significance  in  Q-machine  experiments. 


Figure  2a  represents  the  same  information  as  in  Fig.  1,  but  now 

plotted  as  a  function  of  (Lc/Ls).  Now  the  transition  from  the  local  to  the 

non-local  result  is  seen  to  be  universal,  almost  independent  of  the  shear 

-2 

parameter.  Of  course,  shear  values  with  (p^/Lg)  £  10  will  lead  to  lower 

non-local  values  with  corresponding,  slight,  departures  from  the  universal 

-2 

transition  curve.  For  small  (L  /L  )  <  10  ,  the  case  of  very  small  shear, 

-6  C  S 

p  / L  3  10  °,  retains  the  local  growth  rate  for  the  entire  range 

^  s  -3 

considered,  but  the  case  of  larger  shear,  p^/Lg  =  10  ,  falls  in  the  domain 

L  /p,  <  10  and  the  growth  rate  is  reduced  due  to  filamental  quenching, 

c  1 

Figure  2b  represents  the  same  parameters,  except  for  b  *  1.2  and 

essentially  describes  the  same  phenomena  as  in  Fig.  2a.  We  note  that  both 
the  local  and  the  non-local  limit  growth  rates  are  lower  than  those  for  b  * 


A  typical  solution  for  <p  is  given  in  Figure  3,  for  the  parameters  of 

Fig.  1  in  the  infinite  slab  limit,  L  /p .  +  °°.  The  shear  parameter  is 
-9  ci 

p,/L  =  10  ~  and  the  solutions  obtained  by  the  analytical  method  described 
i  s 

here  or  by  a  direct  numerical  solution  of  the  differential  equation  were 

found  to  be  identical.  The  solution  is  centered  around  u  ~  0.09  and  has 

mo 

a  width  of  iu  ~  0.01.  The  angular  width  decreases  for  lower  shear,  without 
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any  significant  change  in  the  center  of  |$|*  Figures  4a  and  4b 
respectively  represent  the  variations  in  the  complex  angle  around  which 

the  wave  packet  for  the  electrostatic  potential  4>  is  formed,  and  umo  the 
real  angle  where  |$l  has  its  maximum  along  the  real  angle  axis.  The  angle 
unG  experimentally  measurable  in  the  sense  of  normal  mode  dominance  at 
that  angle  as  measured  from  the  local  magnetic  field  orientation.  For  the 
typical  case  of  b  *  0.6,  we  note  the  progression  of  u^  from  0.13  in  the 
local  limit  to  (0.099)  -  i(0.043)  in  the  non-local  (shear  dominant)  limit 
and  a  corresponding  transition  in  uao  from  0.13  to  0.09.  These  curves  are 
drawn  for  p^/L^  33  10  but  they  represent  quite  accurately  the  results  for 
higher  shears,  because  of  the  universality  of  the  transition. 

Figure  5  provides  the  variation  of  y  as  a  function  of  b  for  various  Lc 
and  combinations.  The  top  curve  is  the  local  result  obtained  directly 
from  Q  =  0;  this  corresponds  to  the  case  of  a  uniform  magnetic  field  (zero 
shear).  If  we  introduce  a  small  shear  L^/p^  31  10^  and  a  comparable  current 

width  L  /p .  =  106,  we  essentially  reach  the  shear  dominant  non-local  result 

c  i  - - - 

discussed  in  Ref.  1.  The  transition  is  shown,  with  L  /p  *  10^  and  10^; 

c  i 

the  latter  case  is  already  good  enough  to  almost  approximate  the  local 

result.  Thus  we  make  the  transition  from  the  local  to  the  non-local  result 

-2 

as  Lc/bs  increases  from  10  to  1.  The  lowest  curve  is  for  Lc  «  is  = 
102p^  ang  the  difference  from  the  curve  above  essentially  represents  the 
additional  effect  of  the  increased  shear. 

If  we  use  a  prescribed  rather  than  the  self-consistent  shear,  the 
relation  between  the  angle  and  the  physical  distance  is  simpler, 

u(x)  =  u  +  (x/L  ).  (10) 

o  s 

53  u  +  (L  /L  )£ 
o  c  s 

It  can  be  easily  seen  that  (4b)  reduces  to  (10)  if  £  remains  small,  a 
situation  which  is  valid  when  Lc  >  Ls ,  the  case  when  the  shear  effects 
dominate.  We  have  also  solved  eqs.  (7)  and  (8),  using  eq.  (10),  and  the 
results  do  not  significantly  differ  from  the  case  of  self-consistent  shear. 

We  have  shown  in  this  memo  report  that  the  non-local  effect  of  shear 
is  tempered  by  the  finite  size  of  the  current  slab:  (1)  For  Lc  >  Ls  the 


11 


shearing  of  the  magnetic  field  lines  across  the  length  of  the  slab  is 
sufficient  to  allow  the  non-local,  singular  perturbation  effects  to  take 
place,  and  the  results  of  Ref.  1  are  recovered;  ^2)  For  Lc  <  Ls  on  the 
other  hand,  the  finite  slab  size  is  smaller  than  that  needed  to  provide 
sufficient  angular  space  to  form  a  parabolic  domain  for  Qj,  the  imaginary 
part  of  Q,  of  the  type  of  Fig.  5  of  Ref.  1.  Thus  the  shear  dominated  non- 
local  wave  packet  formation  cannot  occur  unhindered  and  for  Lc  <  10  **La>  * 

the  local  results  are  recovered;  (3)  A  new  effect  occurs  for  <  10  o. , 

1  J 

due  to  the  variation  of  the  current  itself;  thus  filamental  currents  of 
only  a  few  ion  gyroradii  may  not  produce  any  instability. 

The  results  obtained  here  are  quite  significant  for  space-plasma 
applications  where  small  shear  and  finite  width  current  channels  obtain  in 
the  auroral  arc  regions.  Implications  for  various  space  and  laboratory 
plasmas  will  be  described  elsewhere. 
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